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The approximate atomic self-interaction corrections (ASIC) method to density functional theory is 
put to the test by calculating the exchange interaction for a number of prototypical materials, critical 
to local exchange and correlation functionals. ASIC total energy calculations are mapped onto an 
Heisenberg pair-wise interaction and the exchange constants J are compared to those obtained 
with other methods. In general the ASIC scheme drastically improves the bandstructure, which 
for almost all the cases investigated resemble closely available photo-emission data. In contrast the 
results for the exchange parameters are less satisfactory. Although ASIC performs reasonably well 
for systems where the magnetism originates from half-filled bands, it suffers from similar problems 
than those of LDA for other situations. In particular the exchange constants are still overestimated. 
This reflects a subtle interplay between exchange and correlation energy, not captured by the ASIC. 
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I. INTRODUCTION 

Theoretical studies based on density functional the- 
ory (DFT) [TJ 12] have given remarkable insights into 
the electronic and magnetic properties of both molecules 
and solids l3l. In particular, a number of these studies 
attempt to quantitatively describe the magnetic inter- 
action in a broad range of systems including transition 
metals hypothetical atomic chains [HIS], ionic solids 
[Zl [3 |9] , transition metal oxides [lOl |TT] and transition 
metal polynuclear complexes [TH [T3J [TJ] . DFT uses an 
effective single-particle picture where spin symmetry is 
generally broken. For this reason exchange parameters J 
are conventionally extracted by using a mapping proce- 
dure, where total energy calculations arc fitted to a classi- 
cal Heisenberg Hamiltonian [U [TS] . This is then used for 
evaluating the Curie or Neel temperatures, the magnetic 
susceptibility and for interpreting neutron diffraction ex- 
periments 

Notably, the accuracy and reliability of the numerical 
values of the J's depend on the functional used for the ex- 
change and correlation (XC) energy, being this the only 
approximated part of the DFT total energy [17]. Cal- 
culations based on well-known local functionals, namely 
the local density approximation (LDA) and the gener- 
alised gradient approximation (GGA), are successful with 
itinerant magnetism in transition metals [1], but largely 
over-estimates the Heisenberg exchange parameters in 
many other situations [TJ |5J IHl EJ E] . Additional cor- 
rections based on the kinetic energy density (metaCGA) 
[18j marginally improves the agreement with experiments 
[9], although an extensive investigation over several solid 
state systems has not been carried out so far. 

These failures are usually connected to the local char- 
acter of the LDA, which is only weakly modified by con- 
structing XC potentials including the gradient, or higher 
derivative of the charge density. A direct consequence is 
that the charge density is artificially delocalized in space, 
leading to an erroneous alignment of the magnetic bands. 
These are also artificially broadened. A typical example 



is that of NiO , which LDA predicts as Mott-Hubbard in- 
stead of charge-transfer insulator. Clearly a qualitative 
failure in describing the ground state results in an erro- 
neous prediction of the exchange parameters. 

One of the reasons behind the inability of LDA and 
GGA of describing localized charge densities is attributed 
to the presence of the self-interaction error (SIE) [12]. 
This originates from the spurious Coulomb interaction of 
an electron with itself, which is inherent to local func- 
tionals. Hartree-Fock (HF) methods, in the unrestricted 
or spin polarised form, are SIE free and produce sys- 
tematic improved J parameters. However, these meth- 
ods lack of correlation and usually overcorrect. A typ- 
ical example is the antiferromagnetic insulator KNiFs 
for which HF predicts a nearest neighbour J of around 
2 meV [TJ [2ni HH |22J [2S] against an experimental value 
of 8.6 meV [24]. Direct SIE subtraction, convention- 
ally called self-interaction corrected (SIC) LDA, also im- 
proves the results and seems to be less affected by over- 
correction [5| I25j . Similarly hybrid-functionals, which 
mix portions of HF exchange with the local density ap- 
proximation of DFT, perform better than local function- 
als and in several situations return values for J in close 
agreement with experiments [T] |8] . 

It is important to note that both methods based non- 
local exchange or SIC, are computationally demanding 
and thus their application to the solid state remains 
rather limited. It is then crucial to develop practical com- 
putational schemes able to provide a good estimate of the 
exchange parameters for those systems critical to LDA, 
which at the same time are not numerically intensive. 
Based on the idea that most of the SIE originates from 
highly localized states, with a charge distribution resem- 
bling those of free atoms, Vogel et al. [2B] proposed a sim- 
ple SIC scheme where the corrections are approximated 
by a simple on-site term. This method was then gener- 
alized to fractional occupation by Filippetti and Spaldin 
|27j and then implemented in a localized atomic orbital 
code for large scaling by Pemmaraju et al. [28]. Despite 
its simplicity the method has been successfully applied to 
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a number of interesting physical systems including , tran- 
sition metal monoxides [27] [29], silver halides [30], no- 
ble metal oxides [31], ferroelectric materials [?fl 15^ 155] . 
high-k materials diluted magnetic semiconductors 
[55] and also to quantum transport [37] [3S] . 
The method is strictly speaking not variational, in the 
sense that a functional generating the ASIC potential via 
variational principle is not available. However, since typ- 
ically the LDA energy is a good approximation of the 
exact DFT energy, although the LDA potential is rather 
different from the exact KS potential, a "practical" def- 
inition of total energy can be provided. In this work we 
evaluate the ability of this approximated energy in de- 
scribing exchange parameters for a variety of magnetic 
systems. 



where is the LDA Hamiltonian, P.^ip^ir) = 

C(r)(dO and v^'^ = u{[p„];r) + v^;}^^H{plO];r), 
with u and v^i^^^ the Hatree and LDA-XC potential 
respectively. 

In equation Q we have used the fact that at the en- 
ergy minimum the matrix of SIC KS-eigenvalues ef^'^^'~^ 
is diagonalized by the KS orbitals ipn- Importantly such 
minimization scheme can be readily applied to extended 
systems, without loosing the Bloch representation of the 
KS orbitals [431144]- 

The ASIC method consists in taking two drastic ap- 
proximations in equation Q. First we assume that the 
orbitals 0m , that minimize the SIC functional are atomic- 
like orbitals (ASIC orbitals) thus 
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II. THE ATOMIC SIC METHOD 

The seminal work of Perdew and Zunger [19] pioneered 
the modern theory of SIC. The main idea is that of sub- 
tracting directly the spurious SI for each Kohn-Sham 
(KS) orbital Vn- The SIC-LDA [39] XC energy thus 
writes 

occupied 



where E^^-^lp^ p^] is the LDA-XC energy and 6^^^ is the 
sum of the self-Hartree and self-XC energy associated to 
the charge density — IV'nP of the fully occupied KS 
orbital V^JJ 
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^U[p:] + Et,^^[p:,0]. 



(2) 



Here U is the Hartree energy and a is the spin index. 

The search for the energy minimum is not trivial, since 
E^l'^ is not invariant under unitary rotations of the occu- 
pied KS orbitals. As a consequence the KS method be- 
comes either non-orthogonal or size-inconsistent. These 
problems however can be avoided [40] [H] [42] by intro- 
ducing a second set of orbitals 0^ related to the canonical 
KS orbitals by a unitary transformation M 
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The functional can then be minimized by varying both 
the orbitals ip and the unitary transformation A4 , leading 
to a system of equations 
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where ^^^^^^(r) and are the SIC potential and the 
projector associated to the atomic orbital <I>m- Secondly 
we replace the non-local projector with its expecta- 
tion value in such a way that the final ASIC potential 
reads 
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where the orbital occupation (essentially the spin- 

resolved Miilliken orbital population) of <i>„j. 

Note that in the final expression for the potential a 
factor a appears. This is an empirical scaling term that 
accounts for the fact that the ASIC orbital <i> in general 
do not coincide with those that minimize the SIC func- 
tional (jlj. By construction a = 1 in the single particle 
limit, while it vanishes for the homogeneous electron gas. 
Although in general < a < 1, extensive testing [28] 
demonstrates that a value around 1 describes well ionic 
solids and molecules, while a value around 1/2 is enough 
for mid- to wide-gap insulators. In the following we will 
label with ASIC 1/2 and ASICi calculations obtained re- 
spectively with a = 1/2 and a — 1. 

Finally we make a few comments over the total energy. 
As pointed out in the introduction the present theory is 
not variational since the KS potential cannot be related 
to a functional by a variational principle. However, since 
typical LDA energies are more accurate than their corre- 
sponding KS potentials, we use the expression of equation 
([ij as suitable energy. In this case the orbital densities 
entering the SIC are those given by the ASIC orbital <&. 
Moreover, in presenting the data, we will distinguish re- 
sults obtained by using the SIC energy ([T]) from those 
obtained simply from the LDA functional evaluated at 
the ASIC density, i.e. without including the (5„ correc- 
tions 

III. RESULTS 
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^siCpi/) All our results have been obtained with an implemen- 

4^ " " ' tation of the ASIC method YM based on the DFT code 



Siesta |45j . Siesta is an advanced DFT code using pseu- 
dopotentials and an efficient numerical atomic orbital ba- 
sis set. In order to compare the exchange parameters ob- 
tained with different XC functionals we consider the LDA 
parameterization of Cepcrly and Alder [46j, the GGA 
functional obtained by combining Becke exchange [47] 
with Lee- Yang-Parr correlation [?S] (BLYP), the nonem- 
pirical Purdew, Burke and Ernzerhof (PBE) GGA 
and the ASIC scheme as implemented in reference [25] , 

Calculations are performed for different systems crit- 
ical to LDA and GGA, ranging from molecules to ex- 
tended solids. These include hypothetical H-He atomic 
chains, the ionic solid KNiFs and the transition metal 
monoxides MnO and NiO. DFT total energy calculations 
are mapped onto an effective pairwise Heisenberg Hamil- 
tonian 

— ^ ^ JnmSn ' Sm i (9) 

(nm) 

where the sums runs over all the possible pairs of spins. 
In doing this we wish to stress that the mapping is a con- 
venient way of comparing total energies of different mag- 
netic configurations calculated with different function- 
als. In this spirit the controversy around using the spin- 
projected (Heisenberg mapping) or the non-projected 
scheme is immaterial [SI [SHI [5T] . 
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FIG. 1: (Color on-line) H-He-H chains at an inter-atomic dis- 
tance of I.625A. 
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A. H-He chain 

As an example of molecular systems, we consider H- 
He monoatomic chains at a inter-atomic separation of 
1.625 A (see figure [l]). This is an important benchmark 
for DFT since the wave-function is expected to be rather 
localized and therefore to be badly described by local 
XC functionals. In addition the system is simple enough 
to be accessible by accurate quantum chemistry calcula- 
tions. 

As basis set we use two radial functions ( double- C) 
for the s and p angular momenta of both H and He, 
while the density of the real-space grid converges the 
self-consistent calculation at 300 Ry. Here we consider 
all possible Heisenberg parameters. Thus the triangular 
molecule (Fig|l^) has only one nearest neighbour param- 
eter Jfji the 5-atom chain (Fig|l]3) has both first J\.^ and 
second neighbour jj'g parameters, and the 7-atom chain 
(Fig|lJ;) has three parameters describing respectively the 
nearest neighbour interaction with peripheral atoms Jfj, 
the nearest neighbour interaction between the two middle 
atoms J23 ^iid the second neighbour interaction Jfg. 

Following reference [5], accurate calculations based 
on second-order perturbation theory (CASPT2) p2] are 
used as comparison. The quality of each particular func- 
tionals is measured as the relative mean deviation of the 
nearest neighbour exchange parameters only (Jf2; •^12; 



TABLE L Calculated J values (in meV) for the three different 
H-He chains shown in Fig[l] The CASPT2 values are from 
reference [12], while the SIC-B3LYP are from reference [5]. 
The last two rows correspond to J values obtained from the 
LDA energy calculated at the ASIC density. 

"^121 "^23)1 since those are the largest ones 

1 \Ji — jCASPT2| 

^ ^ \^ I jCASPT2| ■ (-'^'^) 

Our calculated J values and their relative b are pre- 
sented in table |T] where we also include results for a fully 
self-consistent SIC calculation over the B3LYP functional 
(SIC-B3LYP) [5.. It comes without big surprise that the 
LDA systematically overestimates all the exchange pa- 
rameters with errors up to a factor 6 for the smaller J 
(J^'g and Jfg) and an average error b for the largest J 
of about 200%. Standard GGA corrections considerably 
improve the description although the J's are still system- 
atically larger than those obtained with CASPT2. Note 
that the results seem rather independent of the particular 
GGA parameterization, with PBE and BLYB producing 
similar exchange constants. This is in good agreement 
with previous calculations [5]. 

SIC in general dramatically improves the LDA and 
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GGA description and our results for ASICi are reason- 
ably close to those obtained with the full self-consistent 
procedure (SIC-B3LYP). This is an interesting result, 
considering that our ASIC starts from a local exchange 
functional, while B3LYP already contains non-local con- 
tributions. We also evaluate the J parameters by using 
the LDA energy evaluated at the ASIC density (last two 
rows in table |l|. In general this procedure gives J's larger 
than those obtained by using the energy of equation ([T]) , 
meaning that the d^'^ contributions reduce the J values. 

It is then clear that the ASIC scheme systematically 
improves the J values as compared to local functionals. 
The agreement however is not as good as the one ob- 
tained by using a fully self-consistent SIC scheme, mean- 
ing that for this molecular system the ASIC orbit als are 
probably still not localized enough. This can alter the 
actual contribution of (5f/^ to the total energy and there- 
fore the exchange parameters. 
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B. Ionic antiferromagnets: KNiFa 

Motivated by the substantial improvement of ASIC 
over LDA, we then investigate its performances for real 
solid-state systems, starting from KNiFs. This is a pro- 
totypical Heisenberg antiferromagnet with strong ionic 
character, a material for which our ASIC approxima- 
tion is expected to work rather well [28j. It is also a 
well studied material, both experimentally [HI [52] and 
theoretically [TJ [HI [HI [32], allowing us extensive com- 
parisons. The KNiFa has cubic perovskite-like structure 
with the nickel atoms at the edges of the cube, flourine 
atoms at the sides and potassium atoms at the center (see 
Figj2]). At low temperature, KNiFs is a type II antiferro- 
magnetic insulator consisting of ferromagnetic (111) Ni 
planes aligned antiparallel to each other. For our calcu- 
lations we use a double-C polarized basis for the s and p 
orbitals of K, Ni and F, a double-^ for the 3c? of K and Ni, 
and a single-C for the 3d of F. Finally, we use 5x5x5 k- 
points in the full Brillouin zone and the real-space mesh 
cutoff is 550 Ry. Note that the configuration used to 
generate the pseudopotential is that of Ni^+, As^3(f . 

We first consider the band-structure as obtained with 
LDA and ASIC. For comparison we also include results 
obtained with LDA+C/ [531 [SI] as implemented in Siesta 
[55j . In this case we correct only the Ni d shell and we 
fix the Hubbard-?/ and Hund's exchange- J parameters by 
fitting the experimental lattice constant (ag = 4.014 A). 
The calculated values are [7=8 eV and J=l eV. The 
bands obtained with the three methods and the corre- 
sponding orbital projected density of states (DOS) are 
presented in figures [3] and [4] respectively. 

All the three functionals describe KNiFa as an insula- 
tor with bandgaps respectively of 1.68 eV (LDA), 3.19 eV 
(ASICi), and 5.0 eV (LDA+U). An experimental value 
for the gap is not available and therefore a comparison 
cannot be made. In the case of LDA and ASIC the gap 
is formed between Ni states, with conductance band bot- 



FIG. 2: (Color on-line) Cubic perovskite structure of KNlFs. 
Color code: blue=Ni, red=F, Green=K. 
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FIG. 3: Band structure for type 11 antiferromagnetic KNiFa 
obtained with a) LDA, b) ASICi and c) LDA+i7 ((7=8 eV 
and J=l eV). The valence band top is aligned at E—Ep—O eV 
(horizontal line). 

tom well described by Cg orbitals. These are progressively 
moved upwards in energy by the SIC, but still occupy the 
gap. Such feature is modified by LDA+U which pushes 
the unoccupied Cg states above the conductance band 
minimum, which is now dominated by K 4s orbitals. 

In more detail the valence band is characterized by a 
low-lying K 3p band and by a mixed Ni-3d/F 2p. While 
the K 3p band is extremely locaUzed and does not present 
substantial additional orbital components the amount of 
mixing and the broadening of the Ni-3(i/F 2p varies with 
the functionals used. In particular both LDA and ASIC 
predict that the Ni 3c? component occupies the high en- 
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FIG. 4: (Color on-line) DOS for type II antiferromagnetic 
KNiFg obtained with a) LDA, b) ASICi and c) LDA+U 
((7=8 eV and J—1 eV). The valence band top is aligned at 
£■=0 eV (vertical line). The experimental UPS spectrum from 
reference [SB] is also presented (thick green line) . The relative 
binding energy is shifted in order to match the K 3p peak. 



ergy part of the band, while the F 2p the lower. For both 
the total bandwidth is rather similar and it is about 9- 
10 eV. In contrast LDA+U offers a picture where the 
Ni-F hybridization spread across the whole bandwidth, 
which is now reduced to less than 7 eV. 

Experimentally, ultraviolet photoemission spec- 
troscopy (UPS) study of the whole KMF3 (M: Mn, Fe, 
Co, Ni, Cu, Zn) series [SB] gives us a spectrum dominated 
by two main peaks: a low K 3p peak and broad band 
mainly attributed to F 2p. These two spectroscopical 
features are separated by a binding energy of about 
10 eV. In addition the 10 eV wide F 2p band has some 
fine structure related to various Ni d multiplets. An 
analysis based on averaging the multiplet structure |56j 
locates the occupied Ni d states at a bounding energy 
about 3 eV smaller than that of the F 2p band. In figure 
|4] we superimpose the experimental UPS spectrum to 
our calculated DOS, with the convention of aligning in 
each case the sharp K 3p peak. 

It is then clear that ASIC provides in general a better 
agreement with the UPS data. In particular both the 
Ni-3(i/F 2p bandwidth and the position of the Fermi en- 
ergy (Ep) with respect to the K 3p peak are correctly 
predicted. This is an improvement over LDA, which de- 
scribes well the Ni-3(i/F 2p band, but positions the K 3p 
states too close to Ep. For this reason, when we align the 
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TABLE II: Calculated J parameters (in meV) and the 
Miilliken magnetic moment for Ni 3d (Pd) in KNiFa. The ex- 
perimental values for J and ao are 8.2 ±0.6 meV and 4.014A 
respectively while the values in brackets are those from refer- 
ence 0. In the table we report values evaluated at the theo- 
retical (Jth and Pd^) and experimental (Jox and P^") lattice 
constant. ASICiy2 and ASICi are obtained from the LDA 
energies evaluated at the ASIC density. 



UPS spectrum at the K 'Sp position, this extends over E-p. 
Finally in the case of LDA+C/, there is a substantial mis- 
alignment between the UPS data and our DOS. LDA+U 
in fact erroneously pushes part of the Ni d mainfold below 
the F 2p DOS, which now forms a rather narrow band. 

We now turn our attention to total energy related 
quantities. In table |ll] we present the theoretical equi- 
librium lattice constant ap and the Heisenberg exchange 
parameter J for all the functionals used. Experimentally 
we have J=8.2± 0.6 meV [H] . The values of ao and J are 
calculated for the type II antiferromagnetic ground state, 
by constructing a supercell along the (111) direction. Im- 
portantly values of J obtained by considering a supercell 
along the (100) direction, i.e. by imposing antiferromag- 
netic alignment between ferromagnetic (100) planes (type 
I antiferromagnet), yield essentially the same result, con- 
firming the fact that the interaction is effectively only 
extending to nearest neighbors. Furthermore we report 
results obtained both at the theoretical equilibrium lat- 
tice constant (Jth) and at the experimental one (Jcx)- 

Also in this case local XC functionals largely overesti- 
mate J, with errors for Jcx going from a factor 8 (LDA) 
to a factor 4.5 (GGA-type). ASIC improves these val- 
ues, although only marginally, and our best agreement 
is found for ASICi, while ASIC1/2 is substantially iden- 
tical to GGA. Interestingly the ASICi performance is 
rather similar, if not better, to that of meta-GGA func- 
tionals [S]. The situation is however worsened when we 
consider J parameters obtained at the theoretical lattice 
constant. The ASIC-calculated ao are essentially identi- 
cal to those from LDA and about 2% shorter than those 
from GGA. Since J depends rather severely on the lattice 
parameter we find that at the theoretical lattice constant 
GGA- functionals perform actually better than our ASIC. 
Finally, also in this case the J's obtained by simply us- 
ing the LDA energies are larger than those calculated by 
including the SIC corrections (see equation [l]) . 
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In general the improvement of the J parameter is cor- 
related to an higher degree of electron localization, in 
particular of the Ni d shell. In table |TT] the magnetic mo- 
ment of the Ni d shell Pd, obtained from the Miilliken 
population, is reported. This increases systematically 
when going from LDA to GGA to ASIC approaching the 
atomic value expected from Ni^+. 

Our best result is obtained with LDA-f-C/, which re- 
turns an exchange of 10.47 mcV for the same U and J 
that fit the experimental lattice constant. This is some- 
how superior performance of LDA-f-L/ with respect to 
ASIC should not be surprising and it is partially related 
to an increased localization. The Ni ions d shell in oc- 
tahedral coordination splits into t2g and eg states, which 
further split according to Hund's rule. The states are 
all filled, while for the Cg only the majority are. By look- 
ing at the LDA DOS one can recognize the occupied t\g 
orbitals (we indicate majority and minority spins respec- 
tively with t and [) at -3 eV, the e| at -2 eV and the 
at about eV, while the empty e| are at between 1 and 
3 eV above the valence band maximum. 

The local Hund's split can be estimated from the e|- 
separation. The ASIC scheme corrects only occupied 
states [53, and therefore it enhances the local exchange 
by only a downshift of the valence band. From the DOS 
of figure |4] it is clear that this is only a small contri- 
bution. In contrast the LDA+C/ scheme also corrects 
empty states, effectively pushing upwards in energy the 
band. The net result is that of a much higher degree 
of localization of the d shell with a consequent reduction 
of the Ni-Ni exchange. This is similar to the situation de- 
scribed by the Hartree-Fock method, which however re- 
turns exchange parameters considerably smaller than the 
experimental value [501 HIl 1121 US] • Interestingly hybrid 
functionals [7j have the right mixture of non-local ex- 
change and electron correlation and produce J's in close 
agreement with the experiments. 

We further investigate the magnetic interaction by 
evaluating J as a function of the lattice constant. Ex- 
perimentally this can be achieved by replacing K with 
Rb and Tl, and indeed de Jongh and Block [58 early 
suggested a d"" power law with a = 12 ± 2. Our calcu- 
lated J as a function of the lattice constant d for LDA, 
GGA, ASICi and LDA-hC/ (f/=8 eV and J=l cV) arc 
presented in figure ??. For all the four functionals in- 
vestigated J varies as a power law, although the calcu- 
lated exponents are rather different: 8.6 for LDA, 9.1 
for GGA, 11.3 for ASICi and 14.4 for LDA-hJ7. This 
further confirms the strong underestimation of the ex- 
change constants from local functionals. Clearly the rel- 
ative difference between the J obtained with different 
functionals becomes less pronounced for small d, where 
the hybridization increases and local functionals perform 
better. Note that only ASICi is compatible with the ex- 
perimental exponent of 12 ± 2, being the one evaluated 
from LDA-|-f7 too large. Importantly we do not expect 
to extrapolate the LDA-I-C/ value at any distance, since 
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FIG. 5: J as a function of the lattice constant for LDA, GGA, 
ASICi and LDA+t/ ((7=8 eV and J=l eV). The symbols are 
our calculate value while the solid lines represent the best 
power-law fit. 

the screening of the parameters U and J changes with 
the lattice constant. 

In conclusion for the critical case of KNiFs the ASIC 
method appears to improve the LDA results. This is es- 
sentially due to the better degree of locahzation achieved 
by the ASIC as compared with standard local function- 
als. However, while the improvement over the bandstruc- 
ture is substantial, it is only marginal for energy-related 
quantities. The main contribution to the total energy in 
the ASIC scheme comes from the LDA functional, which 
is now evaluated at the ASIC density. This is not suf- 
ficient for improving the exchange parameters, which in 
contrast need at least a portion of non-local exchange. 



C. Transition metal monoxides 

Another important test for the ASIC method is that of 
transition metal monoxides. These have been extensively 
studied both experimentally and theoretically and they 
are the prototypical materials for which the LDA appears 
completely inadequate. In this work we consider MnO 
and NiO, which have respectively half-filled and partially- 
filled 3d shells. They both crystallize in the rock-salt 
structure and in the ground state they are both type- 
II antiferromagnetic insulators. The Neel's temperatures 
are 116 K and 525 K respectively for MnO and NiO. In all 
our calculations we consider double-^ polarised basis for 
the s and p shell of Ni, Mn and O, double-^ for the Ni and 
Mn 3d orbitals, and single-^ for the empty 3d of O. We 
sample 6x6x6 fc-points in the full Brillouin zone of both 
the cubic and rhombohedral cell describing respectively 
type I and type II antiferromagnetism. Finally the real- 
space mesh cutoff is 500 Ry. 

The calculated band structures obtained from LDA, 
ASICi/2 and ASICi are shown in figures[6]and[7]for MnO 
and NiO respectively. These have been already discussed 



extensively in the context of the ASIC method [571 [55] 
and here we report only the main features. For both the 
materials LDA already predicts an insulating behavior, 
although the calculated gaps are rather small and the 
nature of the gaps is not what experimentally found. In 
both cases the valence band top has an almost pure d 
component, which suggests these materials to be small 
gap Mott-Hubbard insulators. The ASIC downshifts the 
occupied d bands which now hybridize with the 0-p man- 
ifold. The result is a systematic increase of the band-gap 
which is more pronounced as the parameter a goes from 
1/2 to 1. Importantly, as noted already before [55], the 
experimental band-gap is obtained for a ~ 1/2. 




FT K L rrx K L rr T k l r 



FIG. 6: Calculated band structure for the type 11 anti- 
ferromagnetic MnO obtained from a) LDA, b) ASIC1/2 and 
c) ASICi . The valence band top is aligned at eV (horizontal 
line). 

We then moved to calculating the exchange parame- 
ters. In this case we extend the Heisenberg model to 
second nearest neighbors, by introducing the first (Ji) 
and second (J2) neighbor exchange parameters. These 
are evaluated from total energy calculations for a ferro- 
magnetic and both type II and type I antiferromagnetic 
alignments. Our calculated results, together with a few 
selected data available from the literature are presented 
in table Hill 

Let us first focus our attention to MnO. In this case 
both the J's are rather small and positive (antiferro- 
magnetic coupling is favorite), in agreement with the 
Goodenough-Kanamori rules [59 and the rather low Neel 
temperature. Direct experimental measurements of the 
exchange parameters are not available and the com- 
monly accepted values are those obtained by fitting the 
magnetic susceptibility with semi-empirical methods |60j . 
Importantly this fit gives almost identical first and second 
nearest neighbour exchange constants. In contrast all the 
exchange functionals we have investigated offer a picture 
where J2 is always approximately twice as large as Ji. 
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FT K L rrx K L rrx k l r 



FIG. 7: Calculated band structure for the type 11 anti- 
ferromagnetic NiO obtained from a) LDA, b) ASIC1/2 and 
c) ASICi. The valence band top is aligned at eV (horizon- 
tal line). 



Method 


MnO 




NiO 




Jl J2 Pd 


Jl 


J2 Pd 


LDA 


1.0 2.5 4.49 (4.38) 


13.0 


-94.4 1.41 (1.41) 


PBE 


1.5 2.5 4.55 (4.57) 


7.0 


-86.8 1.50 (1.59) 


ASICi/2 


1.15 2.44 4.72 (4.77) 


6.5 


-67.3 1.72 (1.77) 


ASICi 


0.65 1.81 4.84 (4.86) 


3.8 


-41.8 1.83 (1.84) 


ASICi^2 


1.27 2.65 4.72 (4.77) 


7.1 


-74.6 1.72 (1.77) 


ASICt 


0.69 2.03 4.84 (4.86) 


4.4 


-47.9 1.83 (1.84) 


SEl" 


0.86 0.95 






HF" 


0.22 0.36 






B3LYP= 


0.81 1.71 






PBEO" 


0.96 1.14 






B3LYP'' 




2.4 


-26.7 


HF" 




0.8 


-4.6 


SIC-LDA-^ 




2.3 


-12 


Expt.^ 




1.4 


-19.8 


Expt.a 




1.4 


-17.0 



TABLE 111: Calculated Ji and J2 in meV for MnO and NiO. 
Pd is the magnetic moment of the d shell calculated from the 
type 11 antiferromagnetic phase. Values in bracket are for Pd 
evaluated from the ferromagnetic ground state. ASlCj^j a-nd 
ASIC^ are obtained from the LDA energies evaluated at the 
ASIC density, a) Ref. [60], b) Ref. [6l|, c) Ref. [62], d) Ref. 
[TT], e) Ref. [IS], f) Ref. |M], g) Ref. ^ 



This gives us a reasonably accurate value of Ji for LDA 
and GGA, but J2 is overestimated by approximately a 
factor 2, in agreement with previous calculations [TD]. 

ASIC systematically improves the LDA/GGA descrip- 
tion, by reducing both Ji and J2. This is related to the 
enhanced localization of the Mn d electrons achieved by 
the ASIC, as it can be seen by comparing the Mn d mag- 



netic moments (P^) calculated for different functionals 



(see table III I. Thus ASICi, which provides the largest 



magnetic moment, gives also J's in closer agreement with 
the experimental values, while ASIC1/2 is not very dif- 
ferent from LDA. 

Importantly for half-filling, as in MnO, the ASIC 
scheme for occupied states is fundamentally analogous 
to the LDA+C/ method, with the advantage that the U 
parameter does not need to be evaluated. Finally, at 
variance with KNiFs, it does not seem that a portion of 
exact exchange is strictly necessary in this case. Hartree- 
Fock [HT] results in a dramatic underestimation of the J 
parameters, while B3LYP |62j is essentially very similar 
to ASICi. Curiously the best results available in the lit- 
erature [HI] are obtained with the PBEO functional [55] . 
which combines 25% of exact-exchange with GGA. 

The situation for NiO is rather different. The ex- 
perimentally available data [MJ ES] show antiferromag- 
netic nearest neighbour and ferromagnetic second near- 
est neighbour exchange parameters. The magnitude is 
also rather different with IJ2I > 10 |Ji|. Standard lo- 
cal functionals (LDA and GGA) fail badly and overes- 
timate both the J's by more than a factor 5. ASIC 
in general reduces the exchange constants and drasti- 
cally improves the agreement between theory and exper- 
iments. In particular ASICi gives exchange parameters 
only about twice as large as those measured experimen- 
tally. 

A better understanding can be obtained by looking at 
the orbital-resolved DOS for the Ni d and the O p orbitals 
(figure [8| as calculated from LDA and ASIC. There are 
two main differences between the LDA and the ASIC 
results. First there is an increase of the fundamental 
band-gap from 0.54 eV for LDA to 3.86 eV for ASIC1/2 to 
6.5 eV for ASICi. Secondly there is change in the relative 
energy positioning of the Ni d and O p contributions to 
the valence band. In LDA the top of the valence band is 
Ni d in nature, with the O p dominated part of the DOS 
lying between 4 eV and 8 eV from the valence band top. 
ASIC corrects this feature and for ASIC1/2 the O p and 
Ni d states are well mixed across the whole bandwidth. A 
further increase of the ASIC corrections {a = 1) leads to 
a further downshift of the Ni d band, whose contribution 
becomes largely suppressed close to the valence band-top. 
Thus, increasing the portion of ASIC pushes NiO further 
into the charge transfer regime. 

Interestingly, although ASIC 1/2 gives the best band- 
structure, the J's obtained with ASICi are in better 
agreement with the experiments. This is somehow sim- 
ilar to what observed when hybrid functionals are put 
to the test. Moreira et al. demonstrated [TT] that J's 
in close agreement with experiments can be obtained by 
using 35% Hartree-Fock exchange in LDA. Moreover, in 
close analogy to the ASIC behaviour, as the fraction of 
exact exchange increases from LDA to purely Hartree- 
Fock, the exchange constants decrease while the band- 
gap gets larger. However, while the best J's are obtained 
with 35% exchange, a gap close to the experimental one 




FIG. 8: Calculated orbital resolved DOS for type II anti- 
ferromagnetic NiO obtained with a) LDA, b) ASIC1/2 and c) 
ASICi. The valence band top is aligned at eV. 



is obtained with B3LYP, which in turns overestimates the 
J's. This remarks the subtile interplay between exchange 
and correlations in describing the magnetic interaction 
of this complex material. Finally, it is worth remarking 
that a fully self-consistent SIC [3S] seems to overcorrect 
the J's, while still presenting the erroneous separation 
between the Ni d and O p states. 



IV. CONCLUSIONS 

In conclusions the approximated expression for the 
ASIC total energy is put to the test of calculating ex- 
change parameters for a variety of materials, where local 
and gradient-corrected XC functionals fail rather badly. 
This has produced mixed results. On the one hand, the 
general bandstructure and in particular the valence band, 
is considerably improved and resembles closely data from 
photo-emission. On the other hand, the exchange con- 
stants are close to experiments only for the case when the 
magnetism originates from half- filled shells. For other 
fillings, as in the case of NiO or KNiFs the ASIC im- 
provement over LDA is less satisfactory, suggesting that 
a much better XC functional, incorporating a portion at 
least of exact exchange, is needed. Importantly ASIC 
seems to be affected by the same pathology of hybrid 
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functional, i.e. the amount of ASIC needed for correct- 
ing the J is different from that needed for obtaining a 
good bandstructure. 
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